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POROUS MEDIUM EQUATION TO HELE-SHAW FLOW WITH GENERAL 

INITIAL DENSITY 

INWON KIM AND NORBERT POZAR 


Abstract. In this paper we study the “stiff pressure limit” of the porous medium equation, 
where the initial density is a bounded, integrable function with a sufficient decay at infinity. Our 
particular model, introduced in |PQV| , describes the growth of a tumor zone with a restriction on 
the maximal cell density. In a general context, this extends previous results of Caffarelli-Vazquez 
m and Kim m who restrict the initial data to be the characteristic function of a compact 
set. In the limit a Hele-Shaw type problem is obtained, where the interface motion law reflects 
the acceleration effect of the presence of a positive cell density on the expansion of the maximal 
density (tumor) zone. 


1. Introduction 


In this paper we consider the following degenerate diffusion equation 

pt - V • (pVp) = pG{p) in R” X (0, oo), (1.1) 

where 

p = P^(p):=^p—\ (1.2) 

m — 1 

G' < 0 and G{pm) = 0 for some pm > 0. The model fjl.ip was introduced in [PWI as a model 
problem which describes the growth of cancer cells, with focus on the mechanical aspect of the cell 
density motion. Here the pressure p discourages the overgrowth of the cell density p over some 
critical density pc, which is normalized here as 1. In |PQV| the convergence of the solution p of (11.11) 
and the corresponding pressure variable p was studied in the stiff pressure limff i.e., as m —> oo, in 
the setting of the weak solutions. In the model of a fluid flow, m characterizes the compressibility of 
the fluid with m —>■ oo representing the incompressible limit. It is shown in |PQV| in the setting 
that p and p converges to the limit functions poo and Poo, satisfying the following equations 

- Apoo = G(poo) in H(t) := {poo(-,t) > 0} = {poo{-,t) = 1}, (1.3) 

(Poo)t - V • (pooVpoo) = PooG(poo) in R" X (0, oo). (1.4) 


We mention that, even at a formal level, it is not clear how to derive from (lOll-dTAl) the velocity 
law of the free boundary of the tumor region, 9{poo = !}■ In |PQV| it was conjectured that the 
normal velocity law (11.51) holds for general solutions. This is what we prove, along with the uniform 
convergence of the density variable away from the boundary of the tumor region. Roughly speaking 
we will show the following (see Theorem 11.21 below for the precise statements). 

(a) As m —>■ oo, pm uniformly converges to 1 inside fl{t) and to poe^^^'>* outside of H(t), 

(b) {poo = 1} equals the closure of Ut>o{^{t) x {f}), 

(c) the set Q{t) evolves with the normal boundary velocity (in the viscosity solutions sense) 


_|VPoo|_ 

1 — min[l, 


on 9H(f). 


(1.5) 
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Note that (a) and (b) above imply that < 1 outside of and thus the term min[l, 

in (c) at a boundary point x € dfl{t) is the outer limit of p from the complement of 0(f). Thus (jl.5|) 
coincides with the velocity law conjectured in |PQV| . See Theorem 1 1.2 1 for a more precise statement. 

Note that (c) indicates that p is generically discontinuous across 90(f). Thus proving the con¬ 
vergence result requires keeping track of the pressure variable, which appears to be, at least when 
0(f) has smooth boundary, continuous across 0(f). In terms of the pressure, the equation can be 
written as 

Pt = {m- l)pAp -I- |Vpp + {m- l)pG{p). (1.6) 

Now to state our main result in precise terms, let us denote by pm and Pm the (density and 
pressure) solutions of (11.11) . We will show the convergence of Pm as m —>■ oo to the viscosity solution 
of the following free boundary problem: 


f -Ap = G{p) in {p(-, f) > 0}, 

<1 V = g{x,t)\Vp\ on 9{p(-,f) > 0}, (FB) 

i{p^ = l}c{p(-,f)>0}. 

Here g{x,t) := t)] boundary velocity coefficient, and p^{x,t) := 

is the density in the “exterior” region. We set g = -boo whenever > 1. 

As for the initial data for the free boundary problem (IFBI) . we shall assume that 

Ho C M" open bounded, 9Ho € C^’^, 

p T, ( 1 - 7 ) 

Pq G C(]R") with 0 < Po <1 and pp —>■ 0 as \x\ —>■ oo. 

Note that Pq is the initial density in the “exterior” region, that is, the region outside of Hg. 


Initial data for pm- In terms of the density variable, we would like to show that pm converge to 
p(‘, t) '■= Xn{t) + P^XQ.(tYi where H(f) = {p(-, f) > 0}. To this end we will show that the convergence 
holds locally uniformly for a “well-prepared” initial density po,m approximating the initial density 
function po := +P^Xng- Our approximation is constructed such that the corresponding solution 
Pm is increasing in time (see Lemma HU). As for the general initial data po,m approximating 
Po, the convergence then will hold in norm due to the convergence result for the specific po_m 
iTheorem 11.21) as well as the contraction inequality for pm (14.121) . While we believe that the 
monotonicity of pm is not an essential ingredient of the convergence proof in section 4, it is not clear 
at the moment whether the uniform convergence result obtained in Theorem 11.21 holds for general 
choices of po,™ (see Corollary 14.91) . 


To construct our specific approximation po,m, let us first assume that p^^n satisfies, for some 
d > 0 which is independent of to, 

pI^ g l 1(K") n ci’i(R"), 0 < p^„ < 1 - 5, 

Pom Po locally uniformly as to —>• oo, (1-8) 

to( 1 - (5)'"||T>^p^^||oo0 asm^oo. 

Next suppose 


po.m := max(P„i(po),P^,n) . 

where Pm was introduced in (11.21) . and po is the unique smooth solution of 


-Apo = G(Po) 

Po = 0 


in Ho, 
on R" \ Ho 


(1.9) 
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As we shall see in the lemma below, this will guarantee that pm is monotone increasing in time. 
After we obtain convergence result for this particular approximation of po, we can use contraction 
for solutions of cu to address the case of general po,m- 

Remark 1.1. Given Pq satisfying (11.71) . we can easily define = p^ * Pi/mi where rji/m is the 
standard mollifier with radius 1/m. Such initial data satisfies the assumptions m- Indeed, we 
can easily estimate \\D‘^PQ,m\\oo < ||pf ||li II< Cm^. The rest of (II. 8|) is standard. These 
assumptions, as in |PQV| , are required to prevent the jump singularity of pm over time at < = 0. 

Let us now state the main result in this paper. 

Theorem 1.2. Let p^ solve dm with pQ ra satisfying dEHD-dUl), and let pm be the corresponding 
pressure variable. Then the following holds: 

(a) fTheorem \2. 18\} There is a unique viscosity solution p of (IFBI) with initial data po, where 
Po solves —Apo = G{pq) in flo, and zero otherwise; 

(b) (Lemma \4.4Y b)) {p^ > 1} is contained in the closure of {p > 0}; 

(c) (Corollary o The pressure variable pm locally uniformly converges to p as long as p is 
continuous; 

(d) ( Corollary \4.8\ l Pm locally uniformly converges to p := X{p>Q}+P^X{p=o} away from d {p > 0}. 

(e) (Corollary \2.21\) assuming that p^ is a Lipschitz continuous function, d{p > 0} has zero 
Lebesgue measure in R” x [0,oo). 

(f) (Provosition \5.2\) d\v{-.t] > 0} is of finite perimeter as long as {■, t) < lond{p{-,t) >0}. 

Note that the free boundary motion law in (IFBI) yields (a) a generic discontinuity of p across 
9{p > 0} and (b) a generic discontinuity of p over time when the region {p^ > 1} nucleates. For 
this reason the convergence of pm and Pm as stated appears to be optimal. 

Remark 1.3. Due to the fact that p is nonzero outside of {p > 0}, the set {pm > 0} will degenerate as 
m —>■ oo and will not converge to {p > 0}. But our result iCorollarv 14.81) implies that for any e > 0, 
the set {pm > e} will be a subset of {p > 0} for sufficiently large m. In fact one can characterize 
{p > 0} as 

{p > 0} = {liminfprn > 0}. 

m—¥oo 

As in m we will be using the notion of viscosity solutions, which is based on comparison principle 
with appropriate choices of test functions. In our problem these will be radial functions in local 
neighborhoods with fixed boundaries. In the viscosity solutions theory, this corresponds to the 
usage of second-order polynomials as test functions for nonlinear elliptic equations (see for instance 
[CILj L Therefore the first crucial step in the argument is to prove the above theorem in the radial 
case. When there is no surrounding density, i.e., when Pq = 0, we rely on Barenblatt solutions, a well- 
known family of radially symmetric, compactly supported solutions of the porous medium equation. 
Based on the convergence of these radial solutions we apply the viscosity solution approach to 
obtain the corresponding result in m- On the other hand, when p^ is non-zero, there are no such 
explicit solutions available in the radial setting. The other challenges we face are the possible jump- 
type discontinuity over time of the tumor set {p > 0} due to the free boundary velocity becoming 
infinite in the law dm when the density reaches one, as well as the source term G(p), which each 
prevent the straightforward application of a comparison principle argument between subsolutions 
and supersolutions. 
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Formal derivation of the free boundary motion law. Before we finish this section let us 
present a formal computation indicating the free boundary velocity law (Il 3. Let us write m as 

pt — Ap = pG{p), where p = p^. 

Formally from the definition of p it should be clear that p and the pressure variable converges to 
the same limit Poo as m —>■ oo. Let us also denote the limit density solution as poo, and suppose 
that Poo is discontinuous across V,{t) = {poo{--,t) > 0} = {poo{-,t) = 1}. Again if we take the time 
derivative of the total mass at the formal level, denoting poo = P, Poo = P and and p~ as poo 
inside and outside of n(t), then we have 


pGip) = Pdx = - 


dt 


dt 




p dx ■ 


f {p+)tdx+ f F(p+- p“) + f {p~)tdx 

dn(t) JdQit) dR"-n(i) 

[ Ap+ [ V{p+- p-)dS + [ pG{p) 

Jn{t) Jdn{t) J 

[ [-\Dp\ + V{p+ - p-)]dS + [ pG{p). 

JdQ(t) J 


p dx 


This computation indicates dEll). 


Outline. In section 2 we will prove the comparison principle and uniqueness for the limiting free 
boundary problem (jFBI) . The main results are Theorem 12.141 and Theorem 12.181 They extend the 
comparison and well-posedness results from [P] for the Hele-Shaw problem with a time-dependent 
free boundary velocity coefficient g. The main challenge is to allow for an infinite coefficient depend¬ 
ing on time. This is handled by a shift in time using the fact that the coefficient is nondecreasing in 
time and possesses a certain regularity. In section 3 we show the convergence in the radially symmet¬ 
ric setting with fixed boundary data. Let us mention that we rely on a compactness argument based 
on integral estimates to derive the convergence of the radial solutions in local neighborhoods. Direct 
derivation of convergence using barriers is an interesting open question at the moment. Our integral 
estimates are modified versions from |PQV| due to the presence of fixed boundaries. In section 4 
we prove the convergence result (Corollary 14.81) based on the comparison principle in section 2 as 
well as the radial convergence result in section 3. Lastly, in section 5 we present an estimate on the 
perimeter of the set {p > 0} based on geometric arguments. 

Remark 1.4. Before completion of this paper we learned that similar results were shown by Mellet, 
Perthame and Quiros |MPQ| following a different approach. Their approach relies on integral esti¬ 
mates, while ours relies on pointwise arguments which yield uniform convergence results. We believe 
that both of our approaches have different merits for applications to different contexts. 

2. Notion of solutions and comparison principle 

2.1. Notation. We will follow the notation from [P]. 

Let E <Z for some d > 1. Then USG{E) and LSG{E) are respectively the sets of all upper 
semi-continuous and lower semi-continuous functions on E. For a locally bounded function u on E 
we define the semi-continuous envelopes u*’^ G USG(R‘^) and G LSC{M.‘^) as 

:= inf {v G USG{R‘‘) : v > u on E} , 
u*,E '■= sup {u G LSG{R^) : V < u on E} . 


( 2 . 1 ) 
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Note that m*’'® : —>■ [—c»,c») and t (—c»,c»] are finite on E. We simply write u* and 

M, if the set E is understood from the context. The envelopes can be also expressed as 


*-.E 


(x) = lim sup {u(y) : y G E, \y — x\ < 5} for a; G i?, = —(—m)*’ 




S^O 


Let us review the shorthand notation for the set of positive values of a given function u : E 
defined on a set E C M" x K., 

V,{u; E) := {(x, t) G E : u{x, t) > 0} , E) := {{x, t) G E : u{x, t) < 0} , 


( 2 . 2 ) 


and the closure il{u;E) := il{u;E). For t G K, the time-slices ilt{u]E), Vlt{u\E) and fl1(u]E) are 
defined in the obvious way, i.e., 

flt{u]E) = {x : (x,t) G n{u;E)^ , etc. 

We shall call the boundary of the positive set in E the free boundary of u and denote it r(M; E), i.e., 

r(M;F;) = {dniu;E))nE. 

If the set E is understood from the context, we shall simply write fl(u), etc. 

For given constant t G M we will often abbreviate 

{t < t} := {(x, t) G M” X K : t < t} , etc. 

2.2. Viscosity solutions. We will consider a general problem for the introduction of the notion of 
viscosity solutions. To be more specific, we will define solutions of the problem 

f E{D^u, Du, u) = 0 in {u > 0} , 

( ut—g\Duf = 0 oni9{M>0}. 

We assume that F is a general elliptic operator F{D^u, Du, u) that satisfies the following: There 
exist constants cq, Ci >0 and 0 < A < A such that 

^Aa(^- -ci\p-q\- Co \z-w\< F{M,p,z) - F{N,q,w) 

- ^aa(-^ - N) + ci\p - q\ + Co \z - w\ 

for all M,N G Syrun, p,q G M", z,w G K, where Pucci extremal operators. This 

guarantees that F has the strong maximum principle and Hopf’s lemma; see 0 . Then we need to 
assume that F„ > 0 and that for some pM > 0 

F(0,0,0)<0 and F(0, 0 ,pm)=0. 

Remark 2.1. In the case of (IFBI) we set F(X,p, u) = — trace AT — G{u). 

For the velocity coefficient g : R” x R —>■ (0, oo] we will assume that 

g is continuous at every point of {g < oo} 


and g{x,t) = lim inf _ g{x,t) for all (x,f) 


(2.3) 


As in the previous papers ElP], we define viscosity solutions in two ways: using barriers and 
using test functions. These two notions will be shown to be equivalent, but each has its advantages 
in certain arguments. We will use the notion using barriers, but we still include the notion via test 
functions to show the relation with the original definition in 0. The main difference from [P] is to 
allow for g = +O 0 . 

Before proceeding with the definition of a viscosity solution, we first recall the definition of 
parabolic neighborhood and strict separation used in [P]: 

Definition 2.1 (Parabolic neighborhood and boundary). 

A nonempty set E G R" x R fs called a parabolic neighborhood if E = U r\ {t < t} for some open 
set U C R" X R and some r G R. We say that E is a parabolic neighborhood of {x,t) G R” x R i/ 
(x,t) G E. Let us define dpE ■.= E\E, the parabolic boundary of E. 
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Now we introduce an important concept in the theory, the notion of strict separation. We shall 
use the version introduced in [P], which differs slightly from the one introduced in [K]1. 

Definition 2.2 (Strict separation). Let E C R" x R fee a parabolic neighborhood, and u,v : E ^ 
be bounded functions on E, and let K <Z E. We say that u and v are strictly separated on K with 
respect to E, and we write u v in K, if 

u*’^ < in K r\ n{u; E). 

Remark 2.2. We do not require non-negative functions above, since taking a semicontinuous envelope 
commutes with taking the positive part and 0 < u*’^ = on n{u; E). 

The following lemma was proved in [P] . 

Lemma 2.3 (c.f. (Pj Lemma 2.14]). Suppose that E is a bounded parabolic neighborhood and u, v 
are locally bounded functions on E. The set 

Qu,v,e := {t : u V in E n {t < t}} (2.4) 

is open and Qu,v,e = (—oo,T) for some T S (— 00 , 00 ]. 

2.2.1. Notion via barriers. We build strict barriers for ()2.2|1 . 

Definition 2.3. Let U C M" x R fee a nonempty open set and let (f G be such that Dcf) ^ 0 

on T{(p-,U). We say that (p is a sub-barrier of (12.2|) in U if there exists a positive constant fe > 0 
such that 

(i) F{D^(p, D(j), (p) < —6 in VL{(j)]U), 

(ii) 4>t — g \D<p'^ < —5 on r(((); U). 

A superbarrier is defined analogously by reversing the inegualities in (i)-(ii) and the sign in front of 
S, and requiring additionally that g < 00 on LE{(j)\U). 

Note that it is enough to consider barriers with finite free boundary velocity since we will explicitly 
require in the definition that the positive set of a viscosity solution that the positive set of the solution 
always contains the set where the free boundary velocity coefficient g is infinite. 

Remark 2.4. The Definition 12.31 does not assume 0 > 0, we can always take the positive part later, 
as needed. This does not seem to play a role in the strict separation in Definition 12.21 

The definition of solutions follows. 

Definition 2.4. We say that a locally bounded, non-negative function u : Q —> [0, 00 ) is a viscosity 
subsolution of (1^ on Q if for every bounded parabolic neighborhood E G Q and every superbarrier 
(j) on U such that u ^e 4> on dpE, we also have u ^e 4> on E. 

Similarly, a locally bounded, non-negative function u : Q ^ [0,oo) is a viscosity supersolution of 
(1^ if {g = 00 } r\ Q G D(m*; Q), and for every bounded parabolic neighborhood E G Q and every 
subbarrier tp on U such that (p ^e u on dpE, we also have p -<e u on E. 

Finally, u is a viscosity solution if it is both a viscosity subsolution and a viscosity supersolution. 

Remark 2.5. Since we require {g = 00 } C D(u) for all v G S{g;Q), we also have to address the 
stability of this. That is, 

{g = 00 } C D( inf v) 
veA 

whenever A G S{g- Q). We need that {g = 00 } = int {g = 00 } for this. Then we use subsolutions of 
the elliptic problem in the interior of the positive phase; they give uniform lower bound. 

Remark 2.6. It is not hard to check that if u (v) is a viscosity sub(super)solution of (12.21) then 

—Au < G(u) and — Av > G(v) in {v(-,t) > 0}. 

Remark 2.7. As is standard in the viscosity theory, it is enough to consider only simple cylinders 
with balls as their base as the parabolic neighborhoods E in Definition 12.41 
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2.2.2. Notion via test functions. Similarly to the previous work in mip], we can give an equivalent 
definition of the notion of viscosity solutions via test functions. In the following definitions, Q is an 
arbitrary nonempty parabolic neighborhood. 

Definition 2.5. We say that a locally bounded, non-negative function u \ Q ^ [0,oo) is a viscosity 
subsolution of (12.21) on Q if 

(i) (continuous expansion) 

n(u; (5) n (5 n {t < r} C n(u; Q) r\ {t < t} U {g = oo} for every r > 0, 

(ii) (maximum principle) 

for any (j) £ such that u* — (f has a local maximum at {x, t) € Q D fl{u; Q) in n{u; Q) (~l 
{t < t}, we have 

(ii-1) ifu*{x,i)>0 then F{D'^(j){x,i),D'^(j){x,t),u*{x,i)) <0, 

(ii-2) if u*{x,i) = 0 then either F{D^(j){x,i),D^(fi{x,i),Q) < 0 or Dcj){x,i) = 0 or (j)t{x,i) — 
g(x,t) \D(j)\^ {x,i) < 0. 

Remark 2.8. The condition (i) in Definition [2^ is necessary to prevent a scenario where a “bubble” 
closes instantly; more precisely, a subsolution cannot become instantly positive on an open set 
surrounded by a positive phase, or cannot fill the whole space instantly, unless the expansion of the 
positive phase happens into the set {g = oo}. 

Definition 2.6. We say that a locally bounded, non-negative function u : Q —^ [0, oo) is a viscosity 
supersolution of (I2.2|l on Q if 

(i) (support) 

(i-1) if i^,T) G 12(11,; Q) then {^,t) G 12(11,; Q) for all (^,2) €Q,t>T. 

(i-2) 

{g = oo} n Q C r2(M,; Q). 

(ii) (maximum principle) 

for any (p G such that — (j> has a local minimum at {x,i) G Q m {2 < 2}, we have 
(ii-1) ifu,f{x,i) >0 then F{D^(p(x,t), D^(j){x,i),Ut{x,i)) >0, 

(ii-2) if u,(i, 2) = 0 then either F{D'^(j){x, 2), D^(j){x, 2), 0) > 0 or D(j){x, 2) = 0 or g{x, 2) < oo 
and (ptix, 2) — g{x, i) (x, i) > 0. 

Remark 2.9. As was noted in [P], assumption Definition 12. 6l i-ll is there only to make our life easier. 

Remark 2.10. The closure in the condition Definition 12. 6l i-2l cannot be removed since f2(M,;Q) is 
a (relatively) open set. If at a given time g becomes +oo on an open set outside of Dt(M,) in the 
previous times, then u, is zero on this set. 

Remark 2.11. As is standard in the theory of viscosity solutions, we can require that the test 
functions p are smooth, even polynomials of at most second order in space and first order in time. 
For (ii-2) we can use only radially symmetric test functions. 

The definition of a viscosity solution follows. 

Definition 2.7. We say that a locally bounded, non-negative function u : Q ^ [0, oo) is a viscosity 
solution of (lO) on Q if it is both viscosity subsolution and viscosity supersolution on Q. 

2.3. Equivalence of notions. We now get a result similar to [0 Proposition 2.13]. 

Proposition 2.12. The definitions of viscosity subsolutions (resp. supersolutions) in Definition \2.5\ 
(resp. \2.6\} and in Definition\2.4\ are equivalent. 
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Proof. The direction from Definition 12.51 follows the proof of [Pj Proposition 2.13]. The only detail 
that we have to check is that the supports of a subsolution and a superbarrier stay ordered at the 
crossing time. Since the continuous expansion of subsolution in Definition I2.5l il is valid only in the 
set {g < oo}, we need to use the fact that a for superbarrier in Definition 12.31 satisfies U) C 
{g < oo}. 

We do not have this issue with supersolutions, so the proof is standard. 

The direction from Definition 12.41 to Definition 12.51 and 12.61 is also standard. The continuous 
expansion Definition [13i;i) can be verified by a comparison with radially symmetric barriers. The 
monotonicity of the support of a supersolution Definition 12. 6 l i-ll. an open set at every time, can be 
shown by a comparison with a stationary subbarrier such as = a{c — |a;|^)+ for appropriate 

constants a, c > 0. □ 

With this proposition, we will from now on use the two notions of subsolutions and supersolutions 
from Definition 12.41 and from Definition 12.51 and 12.61 interchangeably. 

2.4. Viscosity solution classes. 

Definition 2.8. For a given function g and a nonempty parabolic neighborhood Q C K" x K. and g 
satisfying (12.31) we define the following classes of solutions: 

• S{g, Q), the set of all viscosity supersolutions of the Ffele-Shaw problem (12.2|) on Q; 

• S_{g,Q), the set of all viscosity subsolutions of (12.2|) on Q; 

• 5 ( 5 , Q) = S{g, Q) n 5 ( 5 , Q), the set of all viscosity solutions of (12.21) on Q. 

2.5. Basic properties of solutions. A subsolution is a subsolution of the elliptic problem on the 
whole space. 

Proposition 2.13. If u € ^ig,Q) for some g and Q then x i-A u*{x,t) is the standard viscosity 
solution of 

< 0 . 

on {x : {x, t) G Q} for every t GM.. 

Similarly, if u G S{g, Q) for some g and Q, then x 1 —>■ u*(x, t) is the standard viscosity solution of 

F{D‘^'if,DiP,f}) > 0 . 

on Q). 

Proof. The proof is analogous to the proof of KB Lemma 3.3]. □ 

2.6. Comparison principle. 

Theorem 2.14. Let Q be a bounded parabolic neighborhood and let gi and g 2 be two velocity coef¬ 
ficients satisfying (lOl) for which there exists f > 0 such that 

g{x,t) := sup gi < _ inf 52 =: g{x,t) for all {x,t) G Q. (2.5) 

Bf(x,t)nQ Bf.(x,t)nQ ~ 

If u G 5 ( 51 , Q) and v G S{g 2 ,Q) such that u -<q v on dpQ, then u <q v in Q. 

2.7. Proof of comparison principle. We can assume that u G USC{Q) and v G LSC{Q). 

We would like to follow the proof of [B Theorem 2.18]. We will use the assumption (12.5p to 
justify the use of sup- and inf-convolutions. 

The structure of the proof is similar to the previous papers p^lKPllP] . with minor modifications 
to allow for the unbounded velocity coefficient. We first regularize the free boundaries of u and v by 
means of the sup- and inf-convolutions over a set of particular shape to guarantee interior/exterior 
ball property in both space and space-time. The set for inf-convolution is decreasing in time to add 
an additional perturbation, by effectively increasing the free boundary velocity of the supersolution. 
Now, if the comparison fails, the regularized solutions must cross. We first show that due to the 
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continuous expansion of the support of u, and the fact that u and v are sub/supersolutions of the 
elliptic problem, this crossing must happen on the free boundary. At the first contact point, the 
boundaries are locally in space. Moreover, the velocity coefficient gi for the subsolution is 
bounded on the neighborhood of this point. At the regular contact point it is possible to define 
weak normal derivatives of the regularized solutions, which must be ordered by Hopf’s lemma. 
Moreover, we can construct barriers to show that the free boundary velocity law is satisfied with 
these weak normal derivatives. An ordering of the free boundary velocities at the crossing point 
with the additional perturbation above then yields a contradiction. Therefore the comparison holds. 
Let us define the crossing time 

to := sup0„,„;(3, (2.6) 

using the set Qu,v,q defined in (12.41) . We observe that u -<q r; in Q is equivalent to to = oo. 

Let us therefore suppose that to < oo and we will show that this leads to a contradiction. 


2.7.1. Regularization. We shall use the standard sup/inf-convolutions to regularize the free bound¬ 
aries at the contact point. We first introduce the open set t) C R" x R for (x, t) e R" x R and 
r > 0 as 

Sr(a:,t) = {( 2 /,s) : (|y-a;| - r)l -f |s - . 

Note that 'E.r{x,t) C if 2r < f. 

Let T > 0 be such that Q C {t < T}. For given 0 < r < f/2, 0 < i5 < ^ we define 

Z{x,t) = sup u, 

Hr {x,t) 

W{x,t)=_ inf V 

^r — St (^!^) 


for (x, t) € Qr with 

Qr ■■= {{x,t) e Q : Erix^t) c Q] . 

Note that Qr is a parabolic neighborhood. 

The following lemma is standard. 


Lemma 2.15. For all r,S > 0 sufficiently small, Z G USC{Qr), W G LSC{Qr), and 

Z -<Q^ W on dpQr- 


For every {x,t) G Qr there exists {Xu,tu) G 'E.r{x,t) C Q and {xv,ty) G 'E.r-st{x,t) such that 

u{xu,tu) = Z{x,t) and v(xv,ty) = W(x,t). 

Moreover, x Z(x,t) is a suhsolution of the elliptic problem on {x : {x,t) G Qr} and x W(x,t) 
is a supersolution of the elliptic problem on flt{W',Qr). 

The support of Z expands continuously in the sense 

Fl{Z-, Qr) r\ Qr r\ {t < t} C fl(Z; Qr) n {t < r} U {g = oo} for every r > 0. 

Similarly, the support of W is nondecreasing, 

if (C, t) G ; Qr) then {f, t) G ^(W ; Qr) for all {f, t) £ Qr, t > t. 


and 


{g = oo} nQr cn{W-,Qr) (2.7) 

Remark 2.16. We can prove a stronger result that actually Z G 5(g; Qr) and W G S(g; Qr), where 
g and g are sup/inf of g over Er, but we actually never need this. 
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Proof. The semicontinuity and existence of points {xu,tu) and (xv,tv) is standard from semiconti¬ 
nuity of u and v. We can choose r < f 12 and S < ^ sufficiently small so that Z and W are strictly 
ordered on dpQr since u and v are strictly ordered on dpQ. 

To check that x !->■ Z{x,t) and x !->■ W{x,t) are a subsolution and a supersolution of the elliptic 
problem in {x : (x, t) G Qr} and ^2t(lV;Qr), respectively, for every t G M, we just need to recall that 
they are the supremum of subsolutions, respectively the infimum of supersolutions, of the elliptic 
problem due to Proposition [2T3] 

The continuous expansion of Z follows from the continuous expansion of u. Indeed, if (^, r) G 
il{Z-Qr) n Qr and g{^,T) < oo, then < oo on Moreover, there exists {fu,tu) G 

n(u; Q) n r) C Bf{f, t). By the continuous expansion of u, we have 

(?«, tu) G fl{u] Q)r\{t < tu}. 

By the definition of the sup-convolution, we conclude that 

(C, r) € n(Z; Qr) n {t < t}. 

To see that the support of W is nondecreasing, suppose that (^, r) G ; Qr). Then by definition 
’E.r-Srif.^x) C Vt{v]Q). Since w is a supersolution, its support is nondecreasing. Definition I2.6l i-lb 
and therefore Er-st{^,t) C Er-Sri^^t) C fi(v;Q) for all t > t. We conclude that (f,t) G ^2(W;Qr) 
for all t > T. 

Finally, if (C,t) G Qr with g{^,T) = oo, then 32 = oo on Bf{^,T). By Definition I2.6l i-2l 
we have Br{^,T) C D(u;(5). Therefore Bp{^,T) n Qr G Q.{W]Qr) for small p > 0 such that 

Bf—pif^x) D ^r_5T-(^,T). □ 

2.7.2. Contact. Let us define the contact time 

i := sup 0 z.W;Qr < io < 00 

where to was introduced as the crossing time in (12. 6 |) . We will show that this leads to a contradiction. 
Lemma 2.17. Z = W = 0 on ^'~{W;Qr) and Z < W on flf(W;Qr). In particular, Z < W on 

QrC{t< t}. 

Proof. Let us denote 

z{x) := Z{x,i), w{x) := W{x,i) 

for X G D := {x : (x,t) = Qr}- Recall that z and w are a subsolution and a supersolution, respec¬ 
tively, of the elliptic problem by Proposition I2.13| The set V := iItj:{W;Qr) is open, and has an 
exterior ball of radius r/2 at every point of its boundary. By (12.7L 'g<g<ooonD\V. We know 
from the definition of the contact time that Il{Z-,Qr) (1 Qr C {t < t} C Ct{W-,Qr). Let y be such 
that Br/ 2 {y) C V^, we must have z = 0 on Br/ 2 {y) H D by the continuous expansion of the support 
of Z and the monotonicity of the support of W in Lemma 12.151 and (12.71) . z is a subsolution of the 
elliptic problem and therefore z = 0 on 5 ^ 72 ( 2 /) D D. By covering D \ D by such balls, we conclude 
that z = 0 on D \ D. 

Let X G V such that z(x) > w{x). We only need to prove that x G dV, and the conclusion then 
follows. Let U be the connected component of V for which x G U. 

We know that z = 0 on dU ("I D from above, and therefore z < w on dU. liU x {t} (~l dpQr ^ 0, 
then the strong maximum principle for the elliptic problem implies that z < m on [/, a contradiction. 

If [7 X {i} C Qr, we have to give a different argument. Let y G U he a. point of maximum of 
z on U. Clearly z{y) > 0. By the interior ball property, there exists f such that y G Br{f) and 
z = z{y) on Br{f,). Since tp = c for c > pm is a supersolution of the elliptic problem on U, the 
strong maximum principle implies z{y) < pm- In particular, z is a strict subsolution of the elliptic 
problem on Br{^). We therefore cannot have w = z on Br{£,). We conclude that z < w on U hy the 
strong maximum principle. □ 
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We know from Lemma [2T^ that Z W vh QrC\ [t < . Therefore due to Lemma [2 .171 we can 

find 

[x,i) & Tl{Z]Qr)r\n’^{W]Qr). 

Due to Lemma [2.151 there exist points 

G ni9D(u;Q) and (x^, G t) n Q). 

We have Sj.(x, i) C D°(it) and Sr(x„, t„) fl Qr C D(Z). Since Z <W on fl {t < t}, we have 
Er-st{x,t) C ft{v) for {x,t) G Er{Xu,tu) n {t < t} . 

By ordering we have 

[Xu ; D — I (x, t) I . 


2.7.3. Free boundary veloeity. Let mz G [—oo, oo] denote the normal velocity of d:ir{x, i) at {xu,tu), 
which can be expressed as 

tu i 

V r 2 - {tu - t) 

Let us define the set 

E . . 

t<i 

Note that E C ft{v) and {xv,tv) G dE. Let mw denote the normal velocity of the boundary of E at 
{xv,tv)- Since D(u) is nondecreasing, we must have mw > 0. But we can also estimate mz — S > mw 
and therefore 


Txiz — d > mw > 0 . 

We conclude in particular that > t > ty. 


2.7.4. Gradients and velocities. Since {g 2 = oo} C D(u) and [xy, ty) G D'^(u), we must have (x^,, ty) G 
{g 2 < oo}. Since (xu,t„), {xy,tv) G ’E.y{x,t) C Bf(x,t), we can estimate 


gi{xu,tu) < _sup gi < _sup gi < _inf _ g 2 < oo. 

B^{x,i) 

Let u be the unit outer normal to {x : (x,f) G Sr(x„,tu)}. We can define the “weak gradients” 


Z(x — hv, t) 
a := limsup- - - 

/i-iOA n 


/3 


:= liminf 
h —^ 0 -|- 


W(x — hu, i) 
h 


Since x is a regular point of the boundary dU, weak Hopf’s lemma implies a < /3, a < oo and /3 > 0. 

As ty > 0, we have enough space to put a barrier above u as in |KP] in a neighborhood of (x„, ty) 
and prove that 


mz < gi{xu,tu)c( < oo. 


Therefore mw < oo. In particular, ty > t — r + 6t. Therefore we have enough space to put a barrier 
under v as in m in a neighborhood of {xy,tv) and prove that 

oo > mz — 5 > mw > g 2 {xy,tv)fd. 

Note that a subbarrier does not need 32 < oo in the complement of its support Definition 12.31 In 
particular, g 2 {xy,tv) < 00 . Putting all together, we have 

mz < gi{xuGu)a < g2ixy,tv)P < mz - S, 


a contradiction. 
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2.8. Well-posedness of (jFBp . We have the following existence and uniqueness result for (|FBIl . 


Theorem 2.18 (Well-posedness). Suppose that flp C M" is a bounded open set with a boundary. 
Moreover, let € C(R") he a function such that 0 < pg < 1, lini|2,|_),oo = 0. Then there 

exists a unique viscosity solution u of (IFBI) with initial support flp o,nd initial density p^ in the 
sense that u is a viseosity solution of (12.21) in Q = R" x (0,oo) with 

g{x,t) . - ; j , tG(o) TT’ (^■®) 

where g = oo if the denominator is 0, and u satisfies the initial condition as 

{x : u*’^{x,0) > 0 } = {x : it,,Q(x, 0 ) > 0 } = flp- 

The solution is unique in the sense that if u and v are two viscosity solutions of (12.21) with the same 
initial data, then 


u 


*,Q 


= V 


= y*,Q- 


(2.9) 


In the proof of the uniqueness in this theorem we also obtain the following version of the com¬ 
parison principle. 


Theorem 2.19. Let Hq Pq be as in Theorem \ 2.18\ Suppose that u is a viscosity subsolution 
and V is a viscosity supersolution of ( 1 ^ in Q = R” X (0, oo) with g as in (12.8|) . with the initial 
data 


Then 


{x : 0) > 0} C Op C {x : x»_q(x, 0) > 0}. 

U*,Q < Q 


We now proceed with the proof of the well-posedness theorem. Let u and v be two solutions of 
(IFBI) on Q = R” x (0, oo) with the given initial data. We want to prove that they must be equal in 
the sense of ( 1 ^ . 

The basic idea is to perturb one of the solutions to create a strictly ordered pair and then apply 
the comparison principle. To apply Theorem 12.141 for a > 1 and cr > 0 we consider the rescaled 
shifted function 


w{x, t) = v(x, at + a). 

Clearly w € S{g2]Q), where 

g2ix) = ag{x,at + cr). 

We want to show that we can find f > 0 such that the assumptions of the comparison principle 
Theorem 12.141 are satisfied. 

Lemma 2.20. Suppose that g satisfies the assumptions (12.3|) . g is nondecreasing in time, and 
{g = oo} is the epigraph of a function r : R” —> R such that r is continuous at every point in 
{t < oo}. Then for every compact set E C R” x [0,oo), a > 1 and a > 0 there exists r > 0 such 
that ag(x, at + s) > g{y, s) whenever (x, t), {y, s) € E and |(x, t) — {y, s)| < r. 

Proof. Let us set /(x, t) := ag{x, at + a). 

1. We first show that we can find K > 0 such that 


:= dist ({(7 > K} , {/ < oo} (~) FI) > 0. 

Indeed, suppose that (5i = 0 for any k G Af. Thus we can find sequences (xfc,tfc) G {g > k}, 
{yk,Sk) G {/ < oo} (~) E such that \{xk,tk) — {yk,Sk)\ < By compactness of E we can assume 
that up to a subsequence {xk,tk) —>■ {x,t) and {yk,Sk) —>■ (x,t) for some (x,t) G E. In particular, 
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0 < i < c». Since we have g{x,i) > liminik^ao g{xk,tk) = c» by (I2.3p . we deduce T{t) < t. 
Furthermore, as ask + <J < T[yk), continuity of t yields 

at + (7 < t(x) < t, 

a contradiction. Therefore we can choose K > 0 such that Si > 0. 

2. Let S 2 := dist({/ = 00 } ,{/ < K} fl E), and we observe that ^2 > 0 due to (I2.3|) and the 
compactness of E. 

Since the set Q ■= {f < K} D E C {g < K} is compact, we can find a modulus of continuity ui 

of both g on this set, and m := ming min{/, g} > 0. Let us find p > 0 such that U!{p) < (a — l)m. 

We set r := i min{(5i, ^ 2 , p}- 

3. Choose now {x, t), (y, s) G E with |(a:, t) — (y, s)| < r. We now prove that f{x, t) > g{y, s). 

• If f{x,t) = 00 then the conclusion is trivial. 

• If K < f{x,t) < 00 then g{y,s) < K and hence f{x,t) > g{y,s). 

• If K < g{y, s) then f(x, t) = 00 , and again the conclusion is trivial. 

• Finally, if neither of the above is satished, we must have f{x,t) < K and g{x,t) < K. 
Therefore we can estimate using the monotonicity in time and continuity 

f{x, t) = ag{x, at + <j) > ag{x, t) 

= {a- l)g{x, t) + g[x, t) > {a - l)m + g{y, s) - u;{r) 

> 9{y,s). 

This finishes the proof. □ 


Proof of Theorem \2.18[ Uniqueness. Let us first prove uniqueness. Suppose that u and v are two 
viscosity solutions satisfying the inital condition. For simplicity, in the following we write u instead 
of u*’^, and v instead of 11 *,g. 

1. If u is a viscosity solution with initial condition fig, a bounded set, we can compare it with a 
large radially symmetric superbarrier 

Wt = ^M(i?2gl6G(0)i/" _ |2.|2)_ (2.10) 

Indeed, since —>■ 0 as |a;| —>• oo, we can for any T > 0 find R sufficiently large such that 

p{x,t) < i for |a;| > i?, 0 < t < r and fig C Bji{0). Then Wt is a superbarrier for 0 < t < T since 
-AWr = 2G(0) > G(0) > G(Wt) in {Wt > 0}, while 

= 4 > 2 > ^ = g{x, t) on d{WT > 0}. 

The comparison with this superbarrier yields that f}{u;'EP x [0,T]) C S/jeieccoT/n (0) x [0,T]. Let 
us therefore define Qt = i? 2 (Rei 6 G(o)T/n)( 0 ) x [0,r]. 

2. We apply the comparison principle on Qt- Since fig has interior ball condition, by comparison 

with radial subbarriers we can prove that fig CC for t > 0. To see this, consider the function 

w{x, t) = a{{ct + r)^ — \x — a;gp). 


For given 0 < r < 1, we can first choose 0 < a <C 1 such that G(4a) > 2na and then choose 
0 < c <C 1 so that c(c + r)/{2ar‘^) < 1. Then for 0 < t < 1 we have w < Aa and therefore 
G{w) > G(4a) > Aw. Moreover, 


Wt 

\Dw\^ 


2ac{ct + r) c{ct + r) 
Aa^\x — xgp “ 2ar2 


<1<«7 


on d{w > 0}. 


We see that ic is a subbarrier for 0 < t < 1. We conclude that if Og has a interior condition with 
radius r > 0, the free boundary of a solution must expand initially with velocity at least c > 0 given 
above. 

3. Let us fix (T > 0. We can find an open set U C M" with smooth boundary such that 
fig CC G CC flcr{v)- cannot jump outside of fig by the definition of a viscosity solution and 










14 


INWON KIM AND NORBERT POZAR 


therefore ^t(u) C U for alH > 0 sufficiently small. By the strong maximum principle for the elliptic 
problem, we obtain that the solution of the elliptic problem on U with data zero on dU is strictly 
smaller then the solution of the elliptic problem on Since x i-A u{x,t) is a subsolution of the 

elliptic problem on R" for any t > 0, and x i-A v(x, a) is a supersolution of the elliptic problem on 
we conclude that m(-,0) < v{-,a) on 

Let us define w{x,t) = v(x, (1 + a)t + a) for some cr > 0. By the reasoning above, u v on 
dpQx- Lemma [2 . 201 implies that the functions gi = g and g2{x, t) = (1 + cr)g{x, (1 + a)t + cr) satisfy 
the assumptions of Theorem 12.141 on Qx- Therefore u < w on R” x [0,T]. 

Now we send cr —^ 0+ and recover 


< (u*)* < P». 

By shifting u instead of v, that is, considering u{x, (1 + cr)“^(t — cr)) and then sending cr —>■ 0+, we 
also obtain 


u* < {v^y < V*. 

By repeating the same argument with u and v interchanged, we obtain the uniqueness of solutions: 

u* = V*, u* = ?;*. 

Existence. Existence follows from standard Perron-Ishii’s method. We hrst construct appropri¬ 
ate barriers. 

1. Let Zp for p > 0 be the unique solution of the elliptic problem in no+^p(0) with boundary value 

zero, and zero outside of flo + Bp, where Bq{0) = {0}. Since Oq G we see that il,o + Bp{0) G 

for small p > 0 and therefore such Zp exists. Clearly U{x, t) = Zo{x) is viscosity subsolution of (12.2|) 
in R" X (0, oo). 

On the other hand, let us define 

(Zatix) 0 < t < p, 

V{x,t) = \Wi{x,t) p < t < 1, 

\Wk[x,t) k—\<t<k, iteratively = 2, 3,..., 

where Wk (with k = T) was defined in (12.101) . Since p^ < 1 on 90o, g as defined in (12.81) is finite 
in a neighborhood of x {0}. Therefore by continuity, we can find p > 0 sufficiently small and 
a > 0 large enough so that Vr is a viscosity supersolution. 

Note that by continuity of U and V for alH > 0 small, we have 

U*{x,0) = 174a;, 0) = E*(a;,0) = 14(0;, 0) = Zo{x). 

2. Let now u be the supremum of viscosity subsolutions w with initial data w*’‘^(a;, 0) = U{x, 0). 
Since U belongs to this class, we see that u is well-defined and u > U. Moreover, the comparison 
principle, with the perturbation above in the proof of uniqueness, yields 

U*<u*<V\ U^<u<V^. 

In particular, u has the correct initial data. We only need to show that it is a solution. We use 
Definition l2.4l Let us show that u is a subsolution. If not, there exists a parabolic neighborhood and 
a superbarrier which u crosses, even though they are strictly ordered on the parabolic boundary. In 
this case, we can perturb the barrier at the crossing point (making it smaller) and deduce that one 
of the subsolutions must cross the perturbed barrier, leading to a contradiction. 

Similarly, to show that it is a supersolution, we suppose that u crosses a subbarrier, if this 
happens, we can perturb the subbarrier, making it larger, and since the perturbed subbarrier is a 
viscosity subsolution, this makes u larger, contradicting the maximality of u. We therefore only 
need to check that {g = oo} H Q C n(M*_Q; Q). But by our assumption on Pq we have {g = ooj = 
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int {g = oo}. Suppose that Bp{^) x {r} G int^ = oo for some (^, r) and p > 0. Let z be the solution 
of the elliptic problem on Bp{^), and 0 outside of Bp{^). Then 


Z{x, t) 


0 t < r, 
z{x) t > T, 


is a viscosity subsolution. In particular, u, > 0 in Bp{^) x {t > r}. From this we conclude that 
mt \^g — oo^ n CZ and thus {p — oo^ n CZ .q; Q)- 

We have proved that u is the unique solution of (I2.2p with g of the form (12.8p and initial support 
i7o. n 


Corollary 2.21. Suppose Pq is Lipschitz. Then d{p > 0} has Lebesgue measure zero in R" x [0, oo). 


Proof. We will show the following density estimate, which is sufficient to conclude: For any T > 0, 
there exists k = k{T) > 0 such that for any space-time ball Bn+i with radius r(T -|- 2) centered at 
(a;o,to) € d{p > 0} n {t < T}, there is a space-time ball Bn+i of radius kr > 0 which lies in both 
{p > 0} and in Bn+i. 

To show this, let us first prove the ordering 

Pi{x,t) := sup p{x,t) <p 2 {x,t) :=p{x,{l + r)t + r), (2.11) 

\x—y\<kr 

which would hold for an appropriate choice of k if is Lipschitz. To see this, first note that the 
order holds at t = 0, due to the step 2 in the proof of Theorem 12.181 

It is straightforward to check that pi is a viscosity subsolution of (IFBI) with modified normal 
velocity V = \Dpi\gi with 

__1_ 

’ 1 — min[l, pb-I-’ 

where w is the continuity mode of p^, and p 2 is a viscosity supersolution of (iFBl) with normal 
velocity V = \Dp 2 \g 2 , where 

1 + r 

^ 1 — min[l, 

Now suppose Pq is Lipschitz so that uj{s) < Ls for some constant L > 0. Then note that if we 
choose k = then we have gi < g 2 for 0 < r < 1. Therefore the comparison principle 

for (|FBI) yields (14.81) . 

Now to check our original claim, suppose (xq, to) & d{p > 0} (~l {t < T}. Let pi as given in (14.81) . 
then the spatial ball B of radius kr and center Xq lies in the positive set of pi. Due to (14.81) , B also 
lies in the positive set of p at time ti := (1 -I- r)to + r. Due to the monotone increasing nature of 
p, we then end up with a space-time cylinder Bkrixo) x [ti -I- ti -I- kr] lying in the positive set of p. 
Since ti < to + f{T + 1), we can conclude that our density estimate holds. 

□ 


3. Convergence in Local Radial setting 

Here we will introduce the notion of radial solutions and give the convergence proof. To make 
local perturbations of general barriers to make first-order approximations in space and time, we need 
to consider radial barriers with fixed boundaries. 

Definition 3.1. The definition will follow the notion via barriers in Definition \2.4\ which considers 
p outside of the tumor region {p > 0} as given a priori by p^{x,t) = Pq {x)e*^^^^. 

Definition 3.2 (Radial solutions). {(f>, p^) is a radial, classical solution of (IFBI) in the cylindrical 
domain {\x — xq\ < R} x [ti,to] or {|a; — Xq] > R} x [ti,to] if 

P^{x,t) = pf(a;,ti)e(‘"*i)‘^(°\ p^{-,ti) G C^(M") 
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and 

(a) is radial with respect to xq and is smooth in its positive phase; 


(b) (p solves (IFBI) in the classical sense with the free boundary motion law V 

(c) > 0 in \x — xo\ = R for ti < t < to; 


m 


1 - 


.E ’ 




(d) < 1 outside of {(p > 0 }. 

Lemma 3.1. The pair (x{ 0 >o} + X{ 0 =o}P^i'/’) *5 unique pair of functions {p,p) in L°°{Q), 
p e C([0,oo];Li(]R”)), p G Pooip), satisfying 

dtp = Ap + pG{p) inV'{Q), p{0)=p^{0) in 

1 2 ^' 

p = <p on dfl in the sense of trace in W ' (R") for a.e. t > 0 


such that 

( 1 ) p,pGL'-((0,r);Li(R")); 

(2) p{t) is uniformly compactly supported in t € [ 0 ,r]; 

(3) |Vp| G L^Qt); 

(4) dtpeMiQr), dtp&MiQr). 

Here Poo is the Hele-Shaw monotone graph 


Poo (p) 


{ 0 } 0 <p<l, 

[ 0 ,oo) p=l. 


Proof. Let us first prove the uniqueness of the solutions of eu. The statement is analogous to 
|PQV[ Theorem 2.4], with the extra boundary condition for p. 

To apply the Hilbert’s duality method that is outlined in |PQV[ Section 3], we need any two 
solutions (pi,pi) to satisfy 


(Pi - P2)'ipt + (pi 


P 2 )Aip + {piG[pi) - p 2 G{p 2 ))ip dx dt = 0 


(3.2) 


for all ^p G G°°{Qt) with boundary data zero on and at t = T. For ip G G^{Qt) this follows 
m- Then this can be extended to include ip nonzero at t = 0 as in jVj. To extend this to all 
Ip whose support touches the boundary, we need to approximate Aip by Ap, ip G CPP{Qt) in the 
correct norm (at least since pi — P 2 G L°°.) However, this is not possible since Vip yf 0 on the 
boundary in general. We therefore use the fact that Vp G L^(0, T; iL^(H)), and show first 


(Pi - P2)'0t 


V(pi -P 2 ) • Vip + (piG(pi) - p 2 G{p 2 ))ip dx dt = 0, 


by approximation and then integrate the second term by parts in space and use that pi = p 2 on dTl. 
Then we just follow |PQV[ Section 3] since the rest does not see the boundary values. 

To finish the proof, we have to show that {p^,(p) satisfies (13.11) . Let us set p = (p and p = 
X{ 0 >o} +X{</>=o}P 0 • We see that p G P°°{p), (p,p) has all the regularity required by the assumptions 
on {p^,(p), and has the correct initial and boundary data. We therefore only need to show that it 
satisfies iH) in the sense of distributions. Let p G C'^{Qt) be a test function. We will verify that 

/ ppt + pAp + pG[p)p dx dt = 0. 

J Qt 
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Since the boundary 5{p > 0} is assumed to be smooth, its unit outer normal is ~ / 

where V is the normal velocity oi d{p > 0} at the given boundary point. Therefore it follows that 

I P‘ft = - I Ptp- I (1 - p)p-7^== dS, 

Jqt -^{^=0} Ja{p>o} V1 + b 

[ pA<f = ( (fiAp + [ |Vp| p dS, 

JQt J{p>0} Jd{p>0} V1 + 

[ pG{p)p= f G{p)p+ [ pG{0)p, 

J Qt *^{^>0} J {p^O} 

We see that the sum of these terms gives zero. □ 

To avoid an initial layer in the limit m —>■ oo, we need to match the initial data for the m-problems. 
For given radial solution (</>,we therefore define the initial data for prn by first finding <5 > 0 so 
that p® (cc, 0) < 1 — (5 on {4>{-, t) = 0} and then setting 

Po.m(a;) = max [P~^{(t){x, 0)), min(l - 5, p^{x, 0))] . (3.3) 

Such 6 can be found due to the assumption (d) in Definition 13.21 and the continuity of p^. Note 
that with above choice of po,m we have po,m -5^ X{ 0 >o} + X{ 0 =o}Pf in and || Vpo,m||Li < G, and 
Po,m = Pmipo.m) 4> Uniformly. Moreover pm is nondecreasing in time when m is large enough 
according to the following lemma. 

Lemma 3.2. Solution pm of (jl.l|) with initial data po,m given in ()3.3I1 is nondecreasing in time for 
m sufficiently large depending only on ||Ap^||oo. 

Proof. Note that due to the comparison principle, we only need to show that pm is nondecreasing 
at the initial time. By the continuity of p^ we can choose a compact set K G {(j) > 0} such that 
p^ < 1 — 5 on K‘^. We can find mo so that </> > Pm(l — S) on K for m > mg. In particular, when 
m > mo, we see that Pm(po,m(a;)) = 4>ix,0) or po,m{x) = p|(x, 0) <1-5. 

In terms of the pressure variable Pm = PmiPm), Pm satisfies the equation 

^tPm — {vn l)p^(Ap^ T Gfjji)') T |^Pm| • 

We know that </> satisfies Acf + G{(j)) = 0 at t = 0 and therefore dtPm > 0 at the points where 
PmiPO.m) — 

On the other hand, if po,m = P^ <1 — 5 then at such point (x, 0) we have 

Apo.m = rn{m - 2)p™“^|Vpo,mP + mp™“^Apo,m > mp™“^Apo,m. (3.4) 

Since po,m < Lm(l — 5) C 1 we have G(po,m) > for sufhciently large m. Since ||Ap^||oo < G, 
we have from (EH) 

Apo^rn + G{po^rn) > —m(l — 5)”* ^||Apo,m||oo + G(0)/2 > 0 

for sufficiently large m uniformly in x. We conclude that dtPm > 0. Therefore pm is nondecreasing 
for sufficiently large m. □ 

Theorem 3.3. For a given radial solution (f> on {|x — xo| < i?} x (ti, <o) or {|x — xo| > R} x (G, to), 
the corresponding solutions Pm, Pm of (HH on the same domain with initial data pm(‘,ti) = P^ at 
t = 0 and boundary data pm = (j) on \x\ = R satisfy the following: Pm uniformly converges to 4>, and 
Pm uniformly converges to pj away from the support of (p. 
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Proof. We consider the case of exterior domain. We will for simplicity assume that xq = 0, ti = 0 
and to = oo. Let T > 0. Let = {i?< |a;| <7-} for r ^ R large enough so that it contains the 
support of solutions Pm for t < T. We shall assume that (j) is smooth up to the boundary dVl and 
> 0 on dfl X [0,oo). Then ())i > 0 in {(/) > 0}. We set 

K = min {^(x, t) : |x| = i?, t > 0} > 0. 

Let us consider the solution pm, Pm = PmiPm), of the porous medium equation (11.11) on 

QT = nx (0,T), 

with initial data at t = 0 and boundary data pm = 4> for |x| = R. 

Estimates. Since (p solves — A$ = G(<1>) in D {■)/)> 0} for every t, we see that dtPmi', 0) > 0. 
Indeed, recall the initial data from (13.31) . For x such that we have 

dtPm = |Vpm|^ >0 at t = 0. 

On the other hand, if Pm{x) = p^{x, 0), we conclude that dtPm > 0 at t = 0 for sufficiently large m 
by the regularity G and the fact that p'f^ < 1. The transition between these two regimes is a 
convex corner (maximum of two nondecreasing initial data). Therefore dtPm > 0 by the comparison 
principle. 

By putting a subsolution under pm, we can hnd i?i /2 > R such that pm(-,t) > k/2 on fti /2 = 
{x : i? < |x| < i?i/2}- 

We first derive the uniform G^’“ and G^’“ estimates for pm on ili/ 2 - Let us rescale in time. Note 
that pm{x,t) := Pm{x, solves the equation 


Pt = pAp + 


m — 1 


|Vpp +pG{p). 


Since p is uniformly away from zero in fli /2 x [0, (m — 1)T] and uniformly bounded from above, 
this is a uniformly parabolic, quasilinear equation in the set considered above. Now uniform G^’“ 
estimate up to the boundary for p, where the G^’“ norm is only depending on the boundary data 
of p as well as the initial data; see Theorem 4.7 and Theorem 5.3 in [Q. We also have uniform G^’“ 
interior estimates up to the initial boundary. In terms of Pm we lose the estimate in time, but we 
still have the estimate in space. Namely, for sufficiently small e > 0 there exists a constant Ct > 0, 
independent of m, such that 


This yields the bound 


+ \\PTni-,t)\\c^’'-{{R+e/2<\x\<R+2e}) < Ct for every 0<t<T. 


\D'^Pm\ + iDpml < G, on {(x,t) : |x| = i? + e, t > 0}. 


(3.5) 


Since the set {x : |x| = i?} is smooth, we can easily create barriers ^ 1 , (j )2 at the boundary that 
coincide with <p on the boundary and pi < p 2 - Moreover, pi is a subsolution and p 2 is a supersolution 
of 


We conclude that 


pt = {m- l)pAp + |Vp|^ + (m - l)pG(p), 


Pi < Pm If p 2 in a neighborhood of {|x| = i?}. 


This will imply that the limit of Pm will have the correct boundary data. 

Uniqueness We shall prove that pm and pm converge to the unique solution of the problem in 
Lemma EH 

The main problem with fixed boundary data arises in the semiconvexity estimate for pm, a variant 
of the Aronson-Benilan estimate. Since the proof relies on the maximum principle for Ap^, we need 
to handle the boundary value of this function. To accomplish this, we use the estimate (13.51) . 
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Indeed, |PQV| derive that w = Apm + G{pm) is a solution of 

wt>{m- l)p,nAw + 2mVpm ■ Via + (m - - (m - 1) {G{pm) - PmG'{pm)) w- (3.6) 

All the arguments here can be made rigorous as explained in lYl Section 9.3]. Since 
minpg[o,pj,f] (G(p) -pG'{p)) > 0, W{t) = - is a subsolution of ([3^. 

Since on F = {(a:, t) : |a:| = i? + e, t > 0} we have (13.51) . we get 


w = Apm + G{pm) > Apm > -G on V 


(3.7) 


for some constant C > 0, independent of m. Let T = sup{t > 0 : IF(t) < —G} = Thus 

W{t) is a subsolution of (|3.6I) with boundary data w{x,t) > W{t) on F (~l {t < T} and therefore 
W (<) < w{x, t) on {0 < t < T}. By a bootstrap argument with a shift W{t — T) for arbitrary r > 0, 
we can deduce that w{x,t) > —G on {{x,t) :\x\> e, t> T}. 

With (13.71) . we can recover all the uniform local L^-estimates on dtPm, ^ Pm, dtPm, ^Pm from 
section 2 of wm , including the L^-continuity of Pm{t) at t = 0. A standard argument implies that 
Pm —>■ P*if, and Pm^ 4’ in x [0,oo)) by the uniqueness result fLemma l3.1D . 

Lipschitz estimate. The functions Pm and pm depend only on r = |a;| and t. In spherical 
coordinates, (EIZl) reads 


Prr T 



Pr + G(j>) > min 




n—1 „ 


We observe that Pr 
7 —Ayr < G we have for Gi = G 

(m—l)i 


.1 Therefore, for given fixed t and all m large so that 

G(0) 


Integration yields 


„l—n 



^Pr) > -Gi. 


r2 ^Pr{r2,t)-r^ V(?’l, 0 > (»'2 “ O > n < T2. 

n 

To get the lower bound on pr{r), r > R e, we use interior parabolic estimates (13.51) which yield 
\pr{R + s,t)\ < G. Therefore 


Pr{r,t) > -G 




r > R + e. 


To get the upper bound, we recall that 0 < p < pm- By the mean value theorem for any r > R 
there exists r 2 G (r, r + I) with \pr{r 2 ,t)\ < pM- Thus 


Pr{r,t)<{^y + 


r + 1 


n —1 


PM + 


Gir 



Therefore pm is locally uniformly Lipschitz in space for every given time t > 0 as long as to > 

G/t + 1. 

Uniform convergence of pm to (p- 

Let us fix Lf C n compact and T > 0. From above we know that Pm 4 in LUm- We can 
find a contable set C {t > 0} dense in {t > 0} and a subsequence of Pm, still denoted by 

Pm, such that Pm{U) 4{4) in L^[K) for every ti. We can choose G = 0 since Pm(uO) —>■ </>(•,0) 
uniformly by the choice of p'p^ in (13.31) . Due to the uniform Lipschitz bound, by taking a subsequence 
if necessary, we can assume that Pm{',ti) —>■ 44,4) uniformly on K for every ti. Let us choose e > 0. 
4 is uniformly continuous on AT x [0,T] and so there exists ^ > 0 such that \4{x,t) — (/)(a;,s)| < e 
for any \t — s\ < 6, x € K. Find TV G N such that Ui^i(^j ~ SjARi + 5/4) D [0,T + 5] and TVf G N 
such that \\pm{-,ti) — 44,ti)\\co < e for alH = 1,..., TV, TO > M. Let now t G [0,r]. We can find 
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1 < i,j < N such that U <t < tj, tj — U < S. Recall that t i-a pm{x,t) is nondecreasing. Thus for 
any x € K and m > M we have 

Pm{x,t) - (j){x,t) < Pm{x,tj) - (l){x,tj) + cj){x,tj) - cj){x,t) < 2e. 

On the other hand 

Pm{x, t) - cj){x, t) > Pm{x, U) - cj){x, t,) + cj){x, t,) - (j>[x, t) > -2e. 

We conclude that the subsequence Pm —t (j) uniformly on K x [0,r]. Since the limit is unique, the 
whole sequence must converge. 

The uniform convergence of pm- We have pm = ■ Let if be a compact 

subset of {p > 0}. But the uniform convergence, there exists e > 0 with pm > e on K for all m 
sufficiently large. Then for every <5 > 0 for all m large we have 

Pm> >1-6. 

The upper bound follows from the uniform upper bound onp^. Therefore pm —t 1 locally uniformly 
in {p > 0}. It remains to show that pm converges to p* locally uniformly away from {p > 0}. 

Lastly we would like to prove the uniform convergence of Pm to p^ outside of {p > 0}. Due to 
the definition of our solution, for each to > 0 there exists ^ > 0 such that < 1 — 5 for some ^ > 0 
outside of {p > 0} for 0 < t < to- Based on this fact we will argue by iteration over small time 
intervals as follows: 

Let us pick T > 0 small and choose xq outside of {p > 0} so that p^{-,T) < 1 — 5 on a 
ring ri < \x\ < 7’2 containing xq. We can pick ri,r 2 such that pm{-,T) uniformly converges to 
p^{-,T) < 1 — (5 at |x| = Pi and |x| = r 2 - This is possible due to the L^-convergence of pm —t p^. 
Since (pm)t > 0, we conclude that Pm stays strictly below 1 — (5 on {|a;| =rj} x [0,r]. 

We now construct a barrier for pm in {ri < |a;| < r 2 } x [0,r] as follows. At t = 0 we pick a 
radial, smooth function <po(x) which has the same value as po near xq and has the value 1 — 5 on the 
boundary |a;| = r^. Now let ip{x,t) = e'^^ipo{x) where C = G(0) + A or something like this. Since 
Pm < 1 — (5 this works fine as a subsolution for the pm equation if T < 0(6), and it follows that pm 
converges uniformly to p^ at Xg- As a consequence we have the uniform convergence of pm to p^ 
in every compact subset of {p > 0}'^ for 0 < T < 0(6). We now iterate over time to conclude up to 
t = tg. Since tg can be chosen arbitrarily large, we conclude. 

□ 

4. Convergence in the general setting 

Based on the Theorem l3.3l next we consider general, i.e., non-radial solutions pm of and the 
corresponding pressure variable Pm = Pm(pm) with initial data po,m given by (II.9|) that approximate 
the initial data dni). 

As we shall see in the lemma below, our choice of initial data po,m will guarantee that pm is mono- 
tonically increasing in time. After we obtain convergence result for this particular approximation of 
Po, we can use contraction for solutions of dm to address the case of general po,m- 

Lemma 4.1. pm increases in time for large enough m. 

Proof. Let us first consider Pm(x,t) := Pg^m(x) exp(tG(0)/2). Writing p = p^ for the sake of brevity, 
we can estimate 

A(p'") + pG(p) = m(m - l)p”^-‘^\Dp\'^ + mp"^-^Ap + pG(p) 

> p (top™“^Ap + G(p)) . 

Due to our assumptions in (11.81) . there exists mg such that the last term is greater than pmG(Q)/2 = 
61tPm and therefore Pm. is a subsolution of dm for m > mg. 
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Additionally, pm{x,t) := P^^{po{x)) is a stationary subsolution of (II.ip . We have defined the 
nondecreasing-in-time functions pm and pm in such a way that max(pm(-, 0), /3m(-, 0)) = /Oo.m- Since 
a maximum of two subsolutions is also a subsolution, we conclude that pm > Txiax(pm, Pm)-, with 
equality at < = 0. Therefore s) > Pmi', 0) for any s > 0. By the comparison principle we have 
Pm(-,s) > Pm(-,t) for any s > t. □ 


Recall that Pm = Pm{p) '■= Our goal is to show their convergence to (IFBI) as m —>■ oo. 

To this end we first define the semi-continuous limits (also referred to as the half-relaxed limits) as 
TO —>■ oo for a family of functions fm as 

liminf* fm{x, t) := lim inf fmix-\-y,t-\- s) 

r->-0 |y| + |s|<r 
m>r~^ 


and 


limsup*:= lim sup fm{x-\-y,t-\-s). 
|y|+hl<’' 

m>r~^ 

Now let us consider the semi-continuous limits of pm and Pm, i-e., 


Pi ■= liminf* pm, Pi '■= liminbpm 


and 


P 2 ■■= lim sup* Pm, P 2 -= lim sup* Pm- 

For technical reasons, it is useful to consider a regularization of p 2 as follows. For a given constant 
(T > 0 let us define 


Pm{x,t):= sup Pm{y,t). 
\y-x\<(j 

Note that p'^ is a subsolution of (HU. Now let us define 


P 2 ■= lim sup* 


Observe that pi is lower semicontinuous and p 2 and pj are upper semicontinuous. Let us also define 
the sets 


Oi(t) := {pi(-,t) > 0}, 02 (t) = {p 2 (-,t) = 1} and = {pj(.,t) = 1}, 

and define p^Cu^) for each t > 0 as the smallest supersolution of —Au = G{u) with Dirichlet 
boundary data in that is, 

P 2 {x,t) '■= inf{w(a:) : w € Arc > G{w) in a domain containing w > 0}, (4.1) 

and we similarly define p 2 corresponding to the set 02 (t). pj is defined in addition to p 2 so that 
we can track the positive set of Pm- P 2 is not sufficient for this purpose since we do not know if pm 
degenerates to zero as to —>• oo inside the set {p 2 = 1}. We use the set ^^(t) instead of 02(t) to 
guarantee that the set is regular enough so that the positive set of P 2 (•, t) coincides with the reference 
set O 2 (t) , as we see in the next lemma. The following lemma shows the relationship between the 
various sets, where the last equality is the only nontrivial relation, and explains the utility of P 2 . 

Lemma 4.2. For any a > 0 we have 

{pi > 0} C {pi = 1} C {pj = 1} = {pj > 0}. 

Proof. Suppose that pi(a;o,to) < 1 for some (a;o,to)- Then there exist mk,Xk,tk, to^ —>■ 00 and 
{xk,tk) (a;o,fo) as k ^ 00 such that Pmk{xk,tk) -?■ pi(a:o,to) < 1- But then pi(xo,to) < 
liminffc^oo jPmk{xk,tk)^ ^ = 0. In particular, {pi > 0} C {pi = 1}. The second inclusion in 
the lemma is due to the fact that pi < P 2 < P 2 for any cr > 0. Lastly, note that due to its definition 
fl 2 (f) is closed and has the interior ball property with balls of radius cr. It now follows from the 
definition of pf that {p^i'jt) > 0} = 


□ 
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We also point out that due to Lemma 14.11 it follows that pi and p 2 are both nondecreasing in 
time. 

Let a;(-) be the mode of continuity for p^. Below we will show that 

(a) p 2 <P 2 (Lemma [43]); 

(b) Pi and P 2 are respectively a supersolution of (IFBI) with p^ = and a subsolution of 

p<j,E ^pE -I- lTheorem l4.6|) : 

(c) pi(-,0) = p2(-,0) is given by (11.31) with Hq ('Lemma 14.71) . 

Due to (b) and the stability property of the viscosity solutions of (IFBI) . we have (P 2 )* < Pi- This 
and (a) yields the convergence results (see Corollary 14.8|) . We first show that (t) (and therefore 
^ 2 {t), fli(t)) is bounded. 

Lemma 4.3. (0 hounded for any t > 0. 

Proof. By our assumption, po uniformly converges to zero as |a;| 00 . Therefore for any T > 0, 

there exists R > 0 such that 

p^{x,0) < for |a;| > R. 

Let us consider the radial solution pm of dnD starting with po,m, where po is given by 

PO = X|x|<R + 

and po,m approximates po as given in Theorem l3.3l Then pm < Pm by comparison principle for (HID. 
Moreover Theorem 13.31 yields that pm uniformly converges to p, which solves (IFBp with p-® < 1/2 
for 0 < t <T, and thus has hnite propagation property up to t = T. Therefore it follows that flj {t) 
is bounded for 0 < t < T and we conclude. □ 

Next we prove the following lemma, to match pfs with p®. 

Lemma 4.4. Let pi,p 2 he as defined above. Then the following holds: 

(a) P 2 < 1 for t >0 and P 2 < M for t > 0; 

(b) Pi > min[l,p®] and {p® > 1} C {pi > 0}; 

(c) P 2 < p”^’^ < 1 outside o/jpf = !}• 


Proof 1. Toshow (a), we write p™ = Pm{pm) andpo.m = Pm{po,m) and we set M := sup^ 2 ,po,m(x), 
which is finite by assumption. By comparison principle, Pm < M for any m. Set 

c = — max G'{s) > 0. 

0<S<M 

The function = Pm + {M — is a supersolution of (11.61) and therefore the 

comparison principle yields Pm < 4>m for all m. P 2 < Tf for t > 0 follows. This then also implies 
Pra < Pm^{4>ra) < max(M,1 for all t > 0. 

2. To show (b), let us choose (xo,to) G K" x (0, oo) and r > 0. To prove the first part, we will 
show that 

pi(a;,t) > min[l,p®(a;o)e®'(°^*] - u}{r)to in Brixo) x [0,to], (4.2) 

where w(r) is the continuity mode of po at xq. Since r > 0 can be arbitrarily chosen and a;(r) —>■ 0 
as r —>■ 0+, we can then conclude. To show (14.21) we consider the function 

<f{x,t) = [a{t)(p{x) — w(r)t]+ for 0 < t < to, 



POROUS MEDIUM EQUATION TO HELE-SHAW FLOW WITH GENERAL INITIAL DENSITY 


23 


where a{t) is an increasing function satisfying a{t) = {po{xo) — w(r)) until it hits and 

a{t) = 1 thereafter, p = tpm is a smoothed version of XBr{xo) that satisfies —Atp"^ < £ for e -C a;(r). 
One can for instance use 

(!;(’■'-I"■ 

We claim that if m is sufficiently large, then (jj satisfies in its positive set 

<4>t- Ap'^ix) 

< G{Q)4> — w(r) + £ 

<G(0)</.-u;(r)/2<Gfy^)</., 

where := ■ Note that the first inequality holds since — > 0 and a{t) < 1, and the 

last inequality holds since (j) < 1 — a;(r). 

Thus (/) is a subsolution of (HID, and it follows from the comparison principle of HU) that (j) < Pm 
and thus (j) < pi, yielding ( 021 ). 

3. Now let us prove the second part of (b) by modifying the subsolution barrier in the above 
step. Suppose p^{-,to) > 1 in Br{xo) for some (a;o,to) and 0 < r < |2G'(0)|“^/^. Since p^ is 
non-decreasing in time, we have p® > 1 on Br{xo) x [to,oo). Then from the first part we have 
Pi > 1 in Br{xo) X [to,oo), and thus for any (5 > 0 and for sufficiently large m{S) we have 

Pm > 1 — (5 in Br{xo) X [tofyi] for > w((5), 

where fy := to + 2G(0)“^(5. 


Now let us construct the barrier (j){x, t) = a(t)ip{x) to compare with pm in Br{xo) x [to,ti], where 

a(t) = e(G'(o)-3i5)(t-ti) 


p{x) 


^ (fy -{x- xo)^) -f (1 - (5) 


so that we have —A((p™) < S and (p > (1 —<5) in Br{xo) with equality on dBr{xo). Also at initial time 
t = to, a{to) = e(G{o)-3S){to-ti) ^ since to — ti = —2G(0)“^(5. Hence we have ip < 1 — 6 < pm 
at t = to and p < 1 — i5 on dBr{xo) x [to,ti]. Also (/) > 1 — 3i5 > ^ in Br{xo) x [tofyi]- 
Then we can estimate 

- A{(^^) <^t- A(^™) 

< [G(0)-3d]<^-hd 

< [G(0) - 5]<k < G{p^^, 

where the first inequality holds due to the fact that a{t) < 1 and — A((p™) > 0 , and the last inequality 
holds for S sufficiently small due to the fact that ^ and p^G'{0) > SG'{0)r^/n > —6/2 for large 
m. Hence we conclude that (j) < pm in Br{xo) x [tofyi]i which yields 

6t^ 

- 3 — <(/>""< Pm in B^/ 2 {xo) X [to,tl] 

8 n ' 

for m > m{S). Thus we conclude that 

Pi{xo,ti) = pi{xo,to + G(0)“U) > 0 (4.3) 


since r is independent of m. As (14.31) holds for arbitrarily small d, it follows that (a;o,to) G {pi > 0} 
and we can conclude. 


4. Lastly to show (c), we will show that for any given (5 > 0 

P 2 < on {P 2 < 1 - 2(5}. 


We will show this iteratively over time intervals of fixed size 7 > 0, where a satisfies 

g(G(0)-Sl)7(l _ J) = 1-5/2. 


(4.4) 

(4.5) 
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Note that (j4.4l) holds for t = 0. Suppose that (j4.4p holds up to t = T, and let us choose (a;o,to) 
in {p 2 < 1 — 25} n {T < t < T + 7 }. Due to the upper-semicontinuity of p 2 and its monotonicity 
in time, there exists r > 0 such that p 2 <^ — 5\Ti i? 2 r(a;o) x Also note that, due to the 

first part of (b) we have min[yO®(•,T), 1] < pi(-,T) < P 2 {-,T) < 1 — <5 < 1 on i? 2 r(a;o) and hence 
^ gG(o)(to-T)p£;(.^ j.) < - (5) < 1 - <5/2 < 1 in B 2 rixo). 

Now based on these facts we will construct a supersolution barrier </ for CH) in E := B 2 r{xo) x 
[T, to] such that </ < in Br{xQ) x [T, to), concluding (I4.4p . 

Let us choose e > 0 and let ip he a. smooth function in B 2 r{xo), (p < I — 5 such that (p = 1 — <5 on 
dB 2 rixQ), p^{-,T) <p< p^{-,T) -|-£ in Br{xo). Now consider the barrier 

4>{x,t) := in E. 

Note that from ()4.5p we have ^<1 — 5/2inE, and thus and thus for large m. Due to 

this fact and that 4> is smooth, it follows that ^ is a supersolution of (HU in E for sufficiently large 
m. Since p 2 <^ — 5hi i? 2 r(a;o) x [T,to], so is for sufficiently large m, and thus p'^ < ^ on 
the parabolic boundary of E. Thus the comparison principle for (HU yields that < </ in E. By 
sending e —>• 0 we conclude that p 2 < p'^’'® at {xo,to), proving (14.4p for the time interval [T,T -I- 7 ]. 
Now we conclude (j4.4ll by iterating our argument over time intervals of length 7 . Lastly we conclude 
(c) by sending <5 —>■ 0 in (14.41) . 

□ 


Next let us prove that p 2 is bigger than the limit supremum of Pm- 
Lemma 4.5. p 2 < P 2 - 

Proof. For any e > 0, take a smooth solution w{x) of —Aw > G{w) + e with 'u; > £ in a domain 
U containing the closure of ^ 2 {to). We will show that p 2 {-,t) < w. Then one can conclude by the 
definition of p 2 ■ 

Due to Lemma [4^ aL (j){x,t) '■= + w{x) is above Pm on the parabolic boundary of 

E := 17 X [ti, to]- 

Moreover, </ is a supersolution of (fra for sufficiently large m since 

+ — l)(j){—A(j) — G{4>)) — > —M/iti — to) + {m— l)e^ — |V?np > 0 for m 1. 

Thus we conclude that pm < </ in E, which yields that p < </ in E. 

□ 


Now we are ready to show our main claim: 

Theorem 4.6. pi and P 2 are respectively a supersolution of (IFBI) with g = a subsolution 

of (|FBP with g’^ = ■ 

First note that Lemma 14.41 will allow us to treat the limiting density outside of the maximal 
density zone essentially as p'®. 

Proof. 1. We will use Definition 12.41 Let us show the subsolution part first. Suppose pj oot 
a subsolution of lEB]) with g'^. This means that there is a superbarrier (f of IFB]) with in 
U ■.= {\x — Zo\ < r} X [ti, 12] which crosses P2 from above at t = to- In other words, we have 

• p 2 (f on the parabolic boundary of 17; 


• p| ^ ^ in 17 (~l {li < 1 < to}; 

• sup|,,_^„|<^(p^ - </)(•, t) > 0 for to<t< t 2 . 
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Since ^ is a superbarrier of (jFBjl . there exists (5 > 0 such that < 1 —2i5 in Bs{xo) x [to — S, to+<5] 
and 

^0 > ^ _ on a{(/)> 0}n{t < to}. (4.6) 

2. From its definition, pj cannot cross ip before its support crosses that of (p. It follows that 
X{p^>o}{-,to) crosses X{ 0 >o} at t = to, and thus along a subsequence > X{ 0 >o}+X{ 0 =o}(p'^’^ + '^) 
for the first time at {xm,tm) with t ti < to as m —>■ oo. Note that the crossing point exist since 
Pm is continuous in time. 

Let xq be a limit point of {xm}- If (pi^Ojti) > 0 then we have a contradiction since in that case it 
can be easily checked that cp is a. supersolution of () 1 . 6 I 1 in a neighborhood of {xo,ti) for sufficiently 
large m. Also due to Lemma 4.4 (c) and the fact that, from Lemma 

{P 2 = 1} = {P 2 (u 0 > 0} C {(pi-,t) > 0} for t < to, 

the limit point {xo,ti) cannot be outside of {(p > 0}. Hence (xo, fi) lies on d{(p > 0}, and ti =to- 

Relying on the continuity of p-®, let us choose 0 < r < i5 such that 

P<t,e ^ p'^^^(^xo,to) + 2 - +5 in D := Br{xo) x [to - r,to + r]. (4.7) 

We now localize (p in H to a radial profile. Since \D(p\ 7 ^ 0 on d{(p > 0}, it follows from the regularity 
of (p that d{(p{-,to) > 0} is a surface. Therefore we can choose r in above definition of D small 
enough such that there is a exterior ball i?r/ 2 ( 2 /o) in {(p{'ito) = 0} touching xo on its boundary. We 
can then use the Taylor expansion of (p up to the second order in space and first order in time to 
construct a new radial superbarrier <p{x,t) = ip{\x — yo|j^) of (IFBI) in D satisfying (14.6|) such that 
{(p(., to) = 0} = B^/ 2 {yo) and tp > (p on the parabolic boundary of D. Then, replacing p with p{x, t) 
for sufficiently small e > 0 if necessary, := Pm{Pm) crosses ip{x,t) in D for large m. Due to 
Lemma I4.4l cl and (14.71) . 

Pv X{^p>o} + [p^’^ ~ 1) > Pm on the parabolic boundary of D (4.8) 

for large m, where p*^’® := p°'’®(a;o, ^ 0 ) + f. 

Now let Pm := Pip,m be the corresponding solution of (11.11) in D, with fixed data 1 on dBr{xo) 
with approximating initial data given as in (13.31) in section [3l Note that, due to the comparison 
principle of dni), Pm ^ Pm in D. On the other hand, the solution (p, p*^’®) of (IFBp in D satisfies 
p ^ ip in D due to (14.6p . Due to Theorem 13.31 limsuD„. - = P'^'^ < P'^'^ + ^ outside of the 

support of p in D, in particular in the zero set of p in D. This contradicts the fact that p^ crosses 
X{v>o} + X{v=o}(p‘^’'® + <^) in D. We can now conclude. 

3. For the supersolution part, first note that the requirement {p® > 1} C {pi > 0} is satisfied by 
Lemma Id^ bl. Next suppose a subbarrier (p of (jFBI) crosses pi from below in {ja; — zo\ > r} x [ti, ^ 2 ] 
at t = to. Parallel arguments as above using Lemma H^ bl would yield the conclusion. 

□ 

Lastly, to apply comparison principle for pi and p 2 , we show that the initial data for p^’s and p^’s 
respectively coincide. 

Lemma 4.7. At t = 0 we have 

(a) limt_,.o+ Pi(-,t) = po := pfxng + XOo locally uniformly away from dXto; 

(b) limt^o+Pi(uO =Po uniformly; 

where po is the unique solution of — Ap = G(p) in Hq with zero boundary data on dflo- 
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Proof. 1. Let us first show (a). First of all note that converges uniformly to po away from 
rio = {po = !}• Also note that, from their definition, converges to 0.2{t) in HausdorfF distance 

as cr —>■ 0. 

Hence by Lemma 14.41 we have 

p^ = Pi = p 2 outside of {p 2 = 1}. (4.9) 

Moreover, by Lemma [4.41 we have pi > 1 on Hq. Hence it is enough to show that 

{P 2 = l}n{t = 0}=n((x{t = 0}. (4.10) 

To this end we consider the domain 

He := {x : d{x, Qq) < 3e} 

for a given e > 0, and choose a point xq £ SHg. By our assumption there exists (5 > 0 depending on 
e such that po < 1 — 25 in i? 2 e(xo), and thus 

p® < 1 — 5 in i? 2 e(a;o) x [0,ti] for some ti > 0. (4-11) 

Let us now consider the radial function 4){x, t) in B 2 e(xo) — B^(t){xo) such that <)> = 0 on dBg(^t'){xo), 
(p = 1 on B 2 eixo) and 

-Ap^x) = G(0) in B 2 eixo) - He(t)(a:o). 

Note that we have \D(j)\ < M/e on dB^t^^^ix^) where M is independent of e as long as e(<) > e/2. 
Combining this fact and (14.111) . it follows that if we choose e(t) = (e— ^t) and p^ (0) = 1 — 25, then 

i4>,p^) is a supersolution of (IFBI) in B 2 eixo) x [Ojfe], where = min[^,fi]. This and Theorem 13.31 
yields that 

P2 < pf < 1 in B^^2{xo) X [0,G]. 

This concludes (|4.10l) and therefore (a). 

2. Next we prove (b). To this end we need to ensure that Pm does not vanish inside of Hq. Again 
this follows from Theorem 13.31 since at each interior point xq £ Qq with Br{xo) £ Hq for some 
r > 0 we can consider radial solution of (IFBI) with p^ = 0 and apply Theorem 13.31 to show that the 
corresponding solutions Pm of (ITel) uniformly converges to ip. Now we can conclude since Pm > Pm 
by the comparison principle of (11.61) . 

3. Now we are ready to prove (b). Fix e > 0 and define 

0/ := {a; : dist(x,R" \ Hq) > e} and Vtg := Hg = {x : dist(x, Hq) < e}. 

In view of (a) and step 2., there exist 5 = 5(e) > 0, to = ^o(e) > 0 and M such that for m> M and 
0 <t < to the following holds: Pm < 5 on dflg , Pm > S in ftf. Let us consider / and g defined by 


and 

Let 


—A/ = G{f) — e in H/ and / = 5 on dflf, 
—Ag = G[g) + e in and p = 5 on 
(p{x,t) := a{t)f{x) and ip{x,t) := b{t)g{x), 


where 

a(t) := min[5e^'^*, 1] and b{t) := max[5“^e“^®'**, 1]. 

Note that the gradient of / is bounded from above in Hy. Using this fact, direct calculations then 
yield that for sufficiently large m and p and ip are respectively subsolution and supersoluton to (II.6|) 
in Hj X (0,to] and x (0,to]- Thus the comparison principle for (11.61) and the choice of 5 and to 
yield 

Ip < Pm in Qg X [0, to] and p^ < </> in H/ x [0, to]. 

Letting m ^ oo and using arbitrarily small e: > 0, we conclude that the Pm’s converge uniformly to 
the solution of the elliptic equation Ho with zero boundary data. 
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□ 


Theorem 14.61 and Lemma 013 together yield our main result: 

Corollary 4.8. Let p be the unique lower-semicontinuous viscosity solution of (IFBI) given by The- 
orem {KTE[ Then the following holds as m ^ oo: 

(a) lim sup* pm = P* and lim inf* pm = P* ■ 

(b) Pm uniformly converges to p := X{p>o} + P^X{p=o} away from d{p > 0}. 

Proof. From Theorem 14.61 and the stability property of viscosity solutions of (|FBI) . it follows that 
p := (liminfo—> 0 ^ 2 )* ^ supersolution of (IFBp with g = Due to Lemma 0771 (a) and the 

convergence of ^^(O) to flo in Hausdorff distance, we conclude that p{-,t) uniformly converges to 
Po{-, 0) as t —>■ 0. 

From the comparison principle Theorem 12.191 it follows that p < pi- Since p 2 < P 2 for any 
(T > 0, it follows that (P 2 )* < P < Pi- Since pi < p 2 by definition, this means (pi)* = (P 2 )* and 
(pi)* = (P 2 )*- This yields that p = pi = (P 2 )* is a viscosity solution of (IFBI) with surrounding 
density p^, and this yields (b). The convergence of pm in the interior of {p > 0} then follows from 
(b). 

It remains to show that pm converges to p^ away from {p > 0}. Note that due to Lemma 1121 

{P > 0} = {P2 > 0} = {p2 = 1}. 

This and Lemma 14.41 (c) yields that limsup*^_,,oo Pm = P 2 < away from {p > 0}. Now we 
conclude by Lemma r4.4f bL which says liminf*^_^ 3 Q pm = Pi > min[l,p^]. 

□ 

Recall that an “almost” contraction property is available for any two solutions pm, Pm of dm 
from |PQV[ (2.12)] in the form 

||pm(i) - Pm(i)||i < e'^(°)‘||pm(0) - Pm(0)||i for any t > 0. (4.12) 

Using the above formula as well as the uniform convergence result obtained in Corollary 14.81 and 
Corollary 2211 we have the following convergence result for general approximating initial data pn^m'- 

Corollary 4.9. Let po := xo,o + PoXO-l 'OJ'ith LIq, Pq as given in (11.71) . with Lipschitz continuous 
Pq . Suppose that po,m converge to po in L^(R"). Then the corresponding solution pm of (II.ip with 
the initial data po.m converges to p as given in Corollary \4-.8\ in the following sense: 

\\Pmit) — p(t)||i —>-0 as m —>• 00 for a.e. t > 0. 

5. A BV ESTIMATE ON THE POSITIVITY SET OE THE PRESSURE 

Here we show that d{p{-, t) > 0} has finite perimeter as long as p^ stays strictly less than 1 near 
d{p{-,f) > 0}. The result already follows from the BV estimates in |PQV| , however our proof is 
based on geometric arguments and thus is of independent interest. 

Lemma 5.1. Let Llt{p) '■= {p(ui) > 0}, where p is as given in Corollary \4.8\ and assume that 
p^ < 1 on dLl(t). Then for given r > 0, there exists sets Llr^t such that 

Llr^t C Llt{p) for each t > 0 

such that 

(a) Llr^t increases with respect to r; 

(b) rir^t has interior ball properties with radius r; 
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(c) |Or,t - Ot(p)| < CeV. 

Proof. To prove this, take the initial positive set 

rig := {x : d{x, fig) > 2r} 

and consider the corresponding approximating solution pm,r of CH) with its limiting initial density 

Po,r ■■= Xiir.o + P^Xqc^- 

Let us now take (^) P 2 defined in (j4.ip with pm,r instead of pm- Let us choose now a = r. 

Then due to (14.6p . P 2 is a subsolution of (IFBI) with g and ^ 2 ( 0 ) = ■ d{x, fig) < r} C flo- Hence 

by comparison principle of (IFBI) we have P 2 ^ 

^2^ = {P2i-^t) > 0} C Ht(p) 

for all t > 0. (c) follows from the contraction inequality 14.121 applied to pm and pm,r given in 

Corollary 14.81 in the limit m —?► 00 . □ 

Proposition 5.2. Under the same assumptions as in Lemma \5.1[ For any r > 0, Qr,t has uniformly 
bounded perimeter. As a consequence {p{-,t) > 0} is a set of finite perimeter. 

Proof. We consider fl" := with r„ = 2“”. We claim that for r < Xn there is at most Cdr^~‘^ 

balls of radius r covering the boundary of ^r„,f 

We will only show the claim for r = Vn, For smaller radius r < r„, the claim holds due to Lemma 
2.5 of [ACM] . We know that increases with respect to n with 

|Or-Or'|<Cr„, (5.1) 

where C is independent of n. Moreover, from the construction above, in fact we have the following 

relation between and 

{a; : d{x, fl'f) < cr^+i} C fl'f+K (5.2) 

where c is independent of the choice of n. 

Now let us take an open covering O of the boundary of 0”“''^ consisting of balls of radius r„+i 
with its center on a boundary point. Let’s take out a family of disjoint balls in O obtained by Vitali’s 
covering Lemma. In each of this disjoint balls, at least one third of the ball is taken by the interior 
of by the interior ball property satisfied at the center of each ball. Also due to (15.21) at least 

a fixed portion of this interior is away from Now we conclude that if the number of the disjoint 
balls are N, then m yields that N{rn+i)‘^ < Cr„+i, or 

N < C(r„+i)i-''. 

Hence we conclude. □ 
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